NMR relaxation rate in non-centrosymmetric superconductors 



K. V. Samokhin 

Department of Physics, Brock University, St. Catharines, 

(Dated: February 2, 2008) 



Ontario, Canada L2S 3A1 



The spin-lattice relaxation rate of nuclear magnetic resonance in a clean superconductor without 
inversion center is calculated for arbitrary pairing symmetry and band structure, in the presence of 
strong spin-orbit coupling. 

PACS numbers: 74.25.Nf, 74.20.Rp, 74.70.Tx 



X 



Most superconducting materials have an inversion cen- 
ter in their crystal lattices. Among a few exceptions 
are the heavy-fermion compounds CePtsSi p] and Ulr 
Q. The peculiar feature of superconductivity without 
inversion center is that there are fewer pairing chan- 
nels available than in a centrosymmetric crystal, because 
the spin degeneracy of the bands is lifted in the pres- 
ence of strong spin-orbit coupling 0, El- Other prop- 
erties of non-centrosymmetric superconductors studied 
recently include paramagnetic limit and Knight shift 
[M Ij l B , I M 0, fToL ITU , helical superconduc ting ph ases 
LL2l LLiJ lidf , and magnetoelectric phenomena [ijj [lj} • 
The determination of the pairing symmetry however 
presents a considerable challenge. Although there are 
strong indications that the gap in CePt3Si has lines of 
nodes, the decisive proof, including the determination of 
the nodal locations, is still lacking. 

A powerful probe of the quasiparticle properties in su- 
perconductors is the spin-lattice relaxation rate T^ 1 of 
nuclear magnetic resonance (NMR). The presence of the 
coherence, or Hebel-Slichter, peak in Tj~ below the criti- 
cal temperature |18| provided an early strong support for 
the Bardeen-Cooper-Schricffer theory of superconductiv- 
ity More recently, the power-law behavior of the 
relaxation rate at low temperatures commonly observed 
in heavy-fermion compounds, see e.g. Ref. [20|, has been 
used as an argument in favor of the existence of gapless 
excitations in those materials. Indeed, if there are line 



see also Ref. |4] , in which the possibility of superconduct- 
ing states with lines of nodes at high-symmetry locations 
appears naturally. In the alternative approach developed 
in Ref. JSj , the spin-orbit coupling is introduced using the 
Rashba model, and the order parameter becomes a mix- 
ture of spin-singlet and spin-triplet components, which 
does not in general have the lines of nodes required by 
symmetry. 

We consider the NMR spin-lattice relaxation due to 
the interaction between the nuclear spin magnetic mo- 
ment "f n I (7„ is the nuclear gyromagnetic ratio) and 
the hyperfme field h created at the nucleus by the con- 
duction electrons. The system Hamiltonian is H = 
H e + H n + Hi n t, where H e describes the electron subsys- 
tem, H n = — 7„ IH is the Zeeman coupling of the nuclear 
spin with the external field if, and Hi nt = —j n Ih is the 
hyperfme interaction. For / = 1/2, we have two nuclear 
spin states I z = ±1/2 with the energies E[ z — —ujoI z , 
where u>q = j n H is the NMR frequency (the spin quan- 
tization axis is parallel to H). Using the Fermi Golden 
Rule, the spin-lattice relaxation rate can be expressed 
[24| in terms of the retarded correlator of the transverse 
components of the hyperfme fields: 
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jgoint) nodes in the gap, then Tf 1 cx T 3 (T 5 ) as T — > T r 

The measurements of the NMR relaxation rate in 
CePtaSi j2^ indicate the likely presence of line nodes. 
On the other hand, a small Hebel-Slichter peak just be- 
low T c observed in Ref. |2S| was not found in other 
experiments [23|. Theoretically, the standard analy- 
sis of the NMR relaxation rate in superconductors, see 
e.g. Ref. 0, is not directly applicable in the non- 
centrosymmetric case because of a complicated spin 
structure of the non-degenerate bands. The purpose of 
the present Rapid Communication is to calculate T^ 1 
in a non-centrosymmetric superconductor with arbitrary 
pairing symmetry. Our analysis does not rely on any spe- 
cific model for the band structure, the only assumption 
being that the bands are well split due to strong spin- 
orbit coupling, which is the case for CePtsSi [3|. In these 
circumstances, it is convenient to use the exact band re 
resentation of the order parameter introduced in Ref. 



where is obtained by an analytic continuation of the 
Fourier transform of the Matsubara correlator (r) = 



h+e " T h-), with h± = h x ±ih y (in our units 

k B = h = i). 

In general, the hyperfme field h can be represented as 
a sum of the Fermi contact, the orbital, and the spin- 
dipolar contributions |24|. Their relative importance de- 
pends on the electronic structure and therefore varies for 
different systems. Let us first assume that the Fermi 
contact interaction is dominant, then h = aS(0), where 
a = — (87r/3)7 e , j e is the electron gyromagnetic ratio, 
and S(r) = (l/2)cr crcr >'4!^(r)^ ( j'(r) is the electron spin 
density (er is a vector composed of Pauli matrices, and 
a =|, I is the spin projection). The relaxation rate Q is 
then proportional to the retarded correlator of the elec- 
tron spin densities at the nuclear site, which can be cal- 
culated by writing the spin operators S± = S x ± iS y in 
the band representation. 

In a non-centrosymmetric crystal with spin-orbit cou- 
pling, the electron bands are non-degenerate almost ev- 
erywhere, except along some high-symmetry lines in the 
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Brillouin zone. Using the Bloch theorem, the spinor field 
operators can be written as 
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where V is the system volume, i is the band index, 
Ci,fc are the destruction operators of band electrons, and 
Xi,k(v , a) — (r,a\i,k) are periodic functions in the unit 
cell, satisfying the normalization condition 
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(f2 is the volume of the unit cell). The band energies 
are ej(fc) = e,(— k). We introduce the notations Ui t k,a- = 

There exist certain relations between the Bloch am- 
plitudes, which follow from the invariance of the sys- 
tem with respect to time reversal. The time-reversal 
operation K = (i<j2)K (K n is the complex conjuga- 
tion) transforms the Bloch spinor corresponding to k in 
band i into a spinor corresponding to — k in the same 
band: K\i,k) = ti(k)\i, —k), where t,(fe) is a non-trivial 
phase factor. From K 2 — — 1 one immediately obtains 
ti(— k) = —ti(k) Q. Then it is straightforward to show 
that 

«i,-fc,T = *i( fc X,fc,i» u i ,-k,l = ~t* (fe)«*, fe , r (3) 

Since the phases of the Bloch functions can be rotated 
independently at different k, there is some freedom in 
choosing Ui^.a and fj(fe). For this reason, no observable 
quantity can depend on ij(fc). 

As an example, let us consider the generalized Rashba 
model, in which the free-electron Hamiltonian is written 
as Hq — eo(fe) + ~f{k)cr, where eo(fe) is an even function 
of k. The effects of spin-orbit coupling are described by 
a real pseudovector "y(k) = — "y(— fc), whose momentum 
dependence is determined by the point symmetry of the 
crystal (the case of the tetragonal group C4v relevant for 
CePt 3 Si is discussed in Refs. HEIII1)- The diagonal- 
ization of Hq gives two Rashba bands with the energies 
e±(fc) = e (fe)±|7(fe)| and M ±, fc , T = ( 1/^2)^1 ± 7z /| 7 |, 
u ±,fe,l = Hlx+ily)/ V 2 l7l(l7l ±Jz)- Applying the time 
reversal operation, one obtains t±(k) = ±.e^^ h \ where 
4> = arg(7 x , j y , 0). It is easy to check that the relations 
(PJ are indeed satisfied. 

The Hamiltonian H e of the electron subsystem consists 
of the free-electron part and the interaction responsible 
for the Cooper pairing. The pairs are formed by electrons 
in the degenerate time- reversed states \i,k) and K\i,k). 
The large band splitting strongly suppresses the pairing 
of electrons from different bands. An anisotropic multi- 
band generalization of the mean-field BCS model reads 

H e = ^ ' &,fc c j,fc c i,fc 
i . k 

+1 X [ A ^ c U4-fc + A? feCi> _ fe c i , fe ] , (4) 



where the chemical potential /i is included in the free- 
electron part, i.e. = £j(fe) — /i. The pairing ampli- 
tudes A^fc = — Ai._fc can be written as = ^(fe)Aj ) fe, 
where the auxiliary gap functions A^fc = Ai._fc all have 
the same symmetry, determined by an even irreducible 
representation T of the crystalline point group. The last 
observation allows us to write ^ = Ylt=i Va 4>a (k), 
where dr is the dimensionality of T, <fra\k) are the even 
basis functions (which can have different fc-dependence 
in different bands while having the same transformation 
properties and sharing the same symmetry-imposed ze- 
ros), and rja are the order parameter components satis- 
fying the mean-field self-consistency equations [25j. 

In the band representation, the spin raising and lower- 
ing operators take the following form: 
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and S-(0) = S\ (0). Inserting these expressions in the 
Matsubara spin correlator one arrives at a two-particle 
Green's function, which can be decoupled in the mean- 
field approximation and represented in terms of the nor- 
mal and anomalous Green's functions for the Hamil- 
tonian (QJ: -(T T Ci t k(r)c\ k ,(0)) = Sij5 k ,k'Gi(k, r) and 

{T T Ci,h(. T ) c i,h'{ty) = SijSk,-k'Fi(k,T). Using the rela- 
tions Q, we arrive at the following expression for the 
spin correlator: 

K+- {V m ) = a2T J2y2 X \ U i,k,l\ 2 \ u 3,k',l\ 2 
n ij,kk f 

x [Gi(k,u n + v m )G 3 {k' ,cu n ) 
+t i (k)t*(k')F}(k,u n + v m )F j (k',u n )], (6) 

where uj n = (2n+l)irT and v m = 2rrnrT are the fermionic 
and bosonic Matsubara frequencies respectively, and 
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FUk,u> n ) = F*(k,uj n ). One can see that the phase fac- 
tors ti(k) drop out of the spin correlator |JBJ, as expected. 

The next steps are standard: first, we perform the 
summation over n, followed by the analytical continua- 
tion K+_(u) = K-\ { v m)\iv m — >w+i0+i then take the ther- 
modynamic limit V — > oo and neglect the electron-hole 
asymmetry, which allows us to write the momentum sums 
as the Fermi-surface integrals. The final result for the re- 
laxation rate QJ is 

oo 

o 

+ Re[Af T »M i H]}, (7) 
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where f(u>) = (e u ' T + 1) 1 is the Fermi function, and 



M (r (cj)=Y J NF } i(\u i ,k,c 
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Here the angular brackets denote the Fermi-surface av- 
erage, and Np,i = (1/87T 3 ) J dSp/\vF,i\ is the density of 
states at the Fermi level in the ith band. The angular 
integration is restricted by the condition (A^kl < lu. 

We see that the intra-band contributions to the NMR 
relaxation rate in a non-centrosymmetric superconductor 
are similar to those in the centrosymmetric case, despite 
a nontrivial spin structure of the single-electron bands. 
The only difference is that it is the auxiliary gap func- 
tions Aj fe that enter Eqs. JSj. It is important however 
that, in addition to the intra-band terms, there are inter- 
band interference contributions to R, which are present 
even without any inter-band scattering mechanisms due 
to electron-electron interactions or impurities 26] . The 
origin of these contributions can be traced back to the lo- 
cal character of the hyperfine coupling IS, which mixes 
together the electron states near the Fermi surface from 
different bands. 

The low-temperature properties of the system are con- 
trolled by the Bogoliubov quasiparticles of energy Ei^ = 

£f k + |Ai.fe| 2 . The experimental data for CePtsSi 

seem to point to the presence of lines of nodes in A^*,. 
Symmetry-imposed gap nodes exist only for the order 
parameters which transform according to one of the non- 
unity representations of the point group. For all such 
non-trivial gap symmetries, the Fermi-surface averages 
of A^fc vanish, and therefore M^(uj) = Mi(oj) = 0. Then 
the relaxation rate in the superconducting normalized to 
its normal-state value can be written as 
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where N n . a — J2i ^F^Ui^^^i are the local densities of 
states of the spin-up (spin-down) electrons at the nuclear 
site r — in the normal state. Note that, although 
|«i,fc,T| 2 + |wi,fc,i| 2 = 1 in the generalized Rashba model, 
this does not have to be the case in general, so that N n ^ + 
N n> i ^ Np, where Np = Yli^F,i is the total density of 
states at the Fermi level for both spin projections. 

Since the order parameters induced by the inter-band 
pair scattering have the same symmetry in all bands, the 
low-energy behavior of all terms in N a (u), see Eq. (JSJ, 
is characterized by the same power law. If there are line 
(po int) nodes in the gap, then N^h o) oc u> (u> 2 ) at u — > 
[111, and R oc T 2 (T 4 ) at T -> [Hlll^. However, a com- 
plicated multi-band structure with considerably different 



local densities of states and gap magnitudes in different 
bands can obscure the simple power-law behavior. In 
particular, in the extreme multi-band case, when some 
of the Fermi-surface sheets remain normal (which might 
be the case in CePtaSi 0), then R = const + aT for line 
nodes, and R = const + aT 2 for point nodes. While the 
residual relaxation comes from the gapless excitations on 
the unpaired sheets, the unusual power laws are due to 
the inter-band contributions to R |2(j. At temperatures 
close to T c , there might still exist the Hebel-Slichter peak 
due to a weak singularity in N a (uj) (for instance, for a 
rf-wave gap the singularity is logarithmic, leading to the 
peak of a much smaller magnitude than in the isotropic 
case) . In general, the presence and the magnitude of the 
Hebel-Slichter peak depend on the gap anisotropy and 
the band structure. 

If the orbital and the spin-dipolar terms are not negli- 
gible, then the calculations become more involved. The 
hyperfine field at the nucleus can be written quite gener- 
ally as a bilinear combination of the band electron oper- 
ators: 



(10) 



ij,kk f 



where A is a pseudovector function, At = A [for the 
Fermi contact interaction considered above, it is factor- 
ized: Aij(k,k') — (a/2)u* k rT er cra rUjk> .&>}■ Undertime 



reversal, c\ k — > U{k)c\_ k and h 
following property: 



-h, which gives the 



Ay(-fe,-fe') = -Ajiik'MUWtffl). (11) 

The symmetry of the Hamiltonian with respect to the 
point group rotations and reflections imposes some addi- 
tional constraints on the A's, which shall not be discussed 
here. 

Assuming that both Ay (fc,fc') and A^fc are weakly 
dependent on ^.fc, using the relations (|ll|l . and repeat- 
ing all the calculation steps from the Fermi-contact case 
above, we obtain 
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(12) 



A^{u) = l\kr 3 {k,k')\ 2 



Re [A* A i)fe ,] 



: , (13) 

AW = (|Ar.( fc) fe')| 2 >... 

The angular brackets here denote the average over the 
Fermi surface in the ith and jth bands, and A ± = 
A^ ± iA v (recall that the spin quantization axis is along 
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the external field H). One can see from Eq. (|f 3fl that 
A^\lo) cx u> 2 (lo 4 ) at lo — > for line (point) nodes, as- 
suming non-vanishing anisotropic gaps in all bands and 
neglecting the possibility of the vertex A^(k, k') going to 
zero accidentally at the gap nodes. Therefore, the low- 
T behavior of the relaxation rate is characterized by the 
familiar exponents: R cx T 2 for line nodes, and R cx T 4 
for point nodes. 

In conclusion, we calculated the NMR relaxation rate 
in a non-centrosymmetric superconductor with strong 
spin-orbit coupling, using the exact band representation 
of the gap functions and the spin density operators. The 
temperature dependence of T-j - is determined by the in- 
terplay of intra-band and inter-band terms, with the lo- 
cal densities of states weighing the contributions from 
different bands. The only robust qualitative conclusion 
that can be drawn from our results is that the low- 
temperature behavior of T^ 1 is determined by the same 
power laws as in the centrosymmetric case: Tf 1 cx T 3 
for line nodes, and T^ 1 cx T 5 for point nodes. 

The determination of the gap symmetry in CePtsSi 
using the NMR relaxation data, beyond the conclusion 
there are likely lines of nodes somewhere on the Fermi 
surface, does not seem to be possible at the moment, 
since it would require a detailed knowledge of the band 



structure and the local densities of states. The absence 
or presence of the Hebel-Slichter peak, especially of the 
magnitude observed in the experiments, cannot be used 
as an argument for or against unconventional supercon- 
ductivity. In addition, the experiments were done at fi- 
nite fields, in the presence of vortices, which might ob- 
scure information about the underlying gap symmetry 
and complicate the comparison with the theory presented 
here. 

Note added: In a very recent theoretical investigation 
of the NMR relaxation rate in CePtaSi [28], a different 
approach was used, in which the spin-orbit band splitting 
is described by the Rashba model, and the gap function 
has both singlet and triplet components. In order to 
explain the observed T 3 behavior of T^ 1 , the parameters 
of the model had to be fine-tuned for the gap to have 
accidental lines of nodes. 



Acknowledgments 

The author is grateful to B. Mitrovic for numerous 
stimulating discussions. This work was supported by 
the Natural Sciences and Engineering Research Council 
(NSERC) of Canada. 



[1] E. Bauer, G. Hilscher, H. Michor, Ch. Paul, E. W. 

Scheidt, A. Gribanov, Yu. Seropegin, H. Noel, M. Sigrist, 

and P. Rogl, Phys. Rev. Lett. 92, 027003 (2004). 
[2] T. Akazawa, H. Hidaka, T. Fujiwara, T. C. Kobayashi, E. 

Yamamoto, Y. Haga, R. Settai, and Y. Onuki, J. Phys.: 

Condens. Matter 16, L29 (2004). 
[3] K. V. Samokhin, E. S. Zijlstra, and S. K. Bose, Phys. 

Rev. B 69, 094514 (2004) [Erratum: 70, 069902(E) 

(2004)]. 

[4] I. A. Sergienko and S. H. Curnoe, Phys. Rev. B 70, 
214510 (2004). 

[5] L. N. Bulaevskii, A. A. Guseinov, and A. I. Rusinov, Zh. 
Eksp. Teor. Fiz. 71, 2356 (1976) [Sov. Phys. - JETP 44, 
1243 (1976)]. 

[6] V. M. Edelstein, Zh. Eksp. Teor. Fiz. 95, 2151 (1989) 

[Sov. Phys. - JETP 68, 1244 (1989)]. 
[7] L. P. Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 

037004 (2001). 

[8] P. A. Frigeri, D. F. Agterberg, A. Koga, and M. Sigrist, 

Phys. Rev. Lett. 92, 097001 (2004). 
[9] P. A. Frigeri, D. F. Agterberg, and M. Sigrist, New J. 
Phys. 6, 115 (2004). 
[10] K. V. Samokhin, Phys. Rev. Lett. 94, 027004 (2005). 
[11] V. P. Mineev, Phys. Rev. B 71, 012509 (2005). 
[12] V. M. Edelstein, J. Phys.: Condens. Matter 8, 339 
(1996). 

[13] K. V. Samokhin, Phys. Rev. B 70, 104521 (2004). 
[14] R. P. Kaur, D. F. Agterberg, and M. Sigrist, preprint 
|cond-m at /0408 149 (unpublished). 



[15] L. S. Levitov, Yu. V. Nazarov, and G. M. Eliashberg, 
Pis'ma Zh. Eksp. Teor. Fiz. 41, 365 (1985) [JETP Lett. 
41, 445 (1985)]. 

[16] V. M. Edelstein, Phys. Rev. Lett. 75, 2004 (1995). 

[17] S. K. Yip, Phys. Rev. B 65, 144508 (2002). 

[18] L. C. Hebel and C. P. Slichter, Phys. Rev. 113, 1504 
(1959). 

[19] M. Tinkham, Introduction to Superconductivity 

(McGraw-Hill, New York, 1996)._ 

[20] J. Flouquet, preprint cond-mat/0501602 (unpublished). 
[21] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
[22] M. Yogi, Y. Kitaoka, S. Hashimoto, T. Yasuda, R. Set- 
tai, T. D. Matsuda, Y. Haga, Y. Onuki, P. Rogl, and E. 

Bauer, Phys. Rev. Lett. 93, 027003 (2004). 
[23] K. Ueda, K. Hamamoto, T. Kohara, G. Motoyama, and 

Y. Oda, preprint cond-mat /0409638 (2004). 
[24] C. P. Slichter, Principles of Magnetic Resonance, 3rd ed. 

(Springer- Verlag, Berlin, 1990). 
[25] V. P. Mineev and K. V. Samokhin, Introduction to 

Unconventional Superconductivity (Gordon and Breach, 

London, 1999). 

[26] K. V. S amokhin and B. Mitrovic, preprint 

cond- mat/0504154 (unpublished). 
[27] If the lines of nodes intersect, then there are logarithmic 

corrections to both N a {uj) and R, see Y. Hasegawa, J. 

Phys. Soc. Jpn. 65, 3131 (1996). 
[28] N. Hayashi, K. Wakabayashi, P. A. Frigeri, and M. 

Sigrist, preprint cond-mat/0504176 (unpublished). 



